Upper critical field of p-wave superconductors with orthorhombic symmetry 
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We extend the Scharnberg-Klemm theory of Hc2 in p-wave superconductors with broken symmetry 
to cases of partially broken symmetry in an orthorhombic crystal, as is appropriate for the more 
exotic ferromagnetic superconductor UCoGe in strong applied magnetic fields. For some partially 
broken symmetry cases, Hc2 can mimic upward curvature in all three crystal axis directions, and 
reasonable good fits to some of the low-field UCoGe data are obtained. 



INRODUCTION 

There has long been an interest in the possibihty of 
superconductivity with the paired electrons having an 
order parameter consisting of a triplet spin configuration 
and the corresponding odd orbital symmetry [1, 2, 3, 4, 
5, 6, 7, 8]. The simplest odd orbital symmetry has the 
p-wave form. In a crystal with a non-cubic structure, 
there can be a variety of different p-wave states [1, 2, 3, 
4, 5] . Depending upon the temperature T, magnetic field 
H, and pressure P, there can be phases corresponding to 
different triplet states [6, 7, 8]. One of the easiest ways 
to characterize the p-wave states is by measurements of 
the T dependence of the upper critical field iJc2 [Ij 2]. 
However, when multiple phases are present in the same 
crystal, as in UPta, a proper analysis requires a variety 
of experimental results [6, 7]. 

Recently, a new class of ferromagnetic superconductors 
has been of great interest. Presently this class consists 
of UGe2[9], Ulr [10], URhGc [11], UCoGe [12], which ex- 
cept for Ulr have orthorhombic crystal structures. For 
URhGe, the superconductivity arises within the ferro- 
magnetic phase. That is also true for UCoGe at ambient 
pressure, but when sufficient pressure is applied, the fer- 
romagnetic phase disappears, leaving the superconduct- 
ing phase without any obvious additional ferromagnetism 
[13, 14]. In the cases of UGe2 and Ulr, applying pressure 
within the ferromagnetic phase induces the superconduc- 
tivity [9, 10]. In addition, polarized neutron experiments 
have been interpreted as providing evidence for a field 
induced ferrimagnetic state in UCoGe, with local mo- 
ments of different magnitudes in opposite directions on 
the U and Co sites [15]. For a ferromagnetic supercon- 
ductor with orthorhombic symmetry, the possible order 
parameter symmetries were given by Mineev [16]. 

Hardy and Huxley measured Bc2{T) of URhGe at am- 
bient pressure in all three crystal axis directions [17]. Us- 
ing only the slope at Tc in each crystal direction as a fit- 
ting parameter, they found that the Scharnberg-Klemm 
theory fit their data quantitatively [17], assuming the 
crystal-aligned polar state with completely broken sym- 
metry (CBS) [2]. For this state, Bc2{T) has a distinctly 
different T dependence for B||c, along the polar axis. 



than for B||d, and B\\b, in the axial plane. By fitting 
only the slopes at in these three crystal axis direc- 
tions, which eliminated the effective mass anisotropy for 
the closed orbits of the paired electrons in the three crys- 
tal directions, as in Landau diamagnetism and in the de 
Haas- van Alphen effect, the results implied that the p- 
wave state oc Pc for all applied field directions. This 
remarkable fit for the low field regime of the supercon- 
ducting state in URhGe did not require the inclusion of 
the ferromagnetism into the theory, as the only apparent 
effect of the magnetization was to give rise to a demag- 
netization jump in Hc2 at the superconducting transition 
temperature, Tc, beyond which the slopes of Hc2{T) were 
fit. In addition, i?c2(0) exceeded the Pauli limit for all 
three crystal axis directions, providing compelling evi- 
dence for parallel-spin state Cooper pairing. 

Upon the discovery of magnetic-field induced reentrant 
superconductivity in URhGe [18], much interest turned 
to the possible source of the high field superconduct- 
ing phase. Then, superconductivity was discovered in 
UCoGe [12], and Hc2{T) was measured for all three crys- 
tal axis directions [19], and all of the curves exhibited 
upward curvature [20]. Subsequently, a highly anoma- 
lous ^'-shaped Hc2{T) curve was observed for T < 0.65Tc 
with H||b. Since M |j c at low fields, this change in 
the M direction only occurred in very pure, well aligned 
samples. This behavior may also have something to do 
with a reentrant phase, one that is close in field strength 
to the low-field phase [21,22]. 

The first attempts to describe upward curvature in 
Hc2{T) in all three crystal directions were based either 
upon ferromagnetic fluctuations [23] , or upon a crossover 
from one parallel-spin state to another [24]. Meanwhile, a 
mean-field theory of the complementary effects of itiner- 
ant ferromagnetism and parallel-spin superconductivity 
was developed [25, 26]. To date, the field dependence 
of this mutual enhancement has not been investigated. 
Here, we study the case in which the p-wave pairing in- 
teraction strength is anisotropic, but finite in all crystal 
directions. Since Hc2{T) is essentially isotropic in the ab 
plane for samples of UCoGe with medium purity [19], we 
studied the partially broken symmetry (PBS) state as a 
function of the pairing interaction anisotropy. This can 
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give rise to a kink in Hc2{T) in at least one crystal axis 
direction [27]. 

CALCULATION 

Our calculation of Hc2{T) assumes weak coupling for 
a homogeneous clean type-II material. Hc2 is therefore 
found by solving the linearized Gor'kov gap equation [2] 

xcxp(— 2^|w„|)exp(— i^mw^sgnw^k ■ M ■ IT) 
xA(R,k'). (1) 
In our present case the interaction is of p-wavc form 

F(k, k') = 3 ^ V^i ■k'ilSi {Xj )><Si{xj)\, (2) 

i 

which can be written in terms of the relevant spherical 
harmonics. Here A''(0) is the single-spin density of states 
at the Fermi energy, 11 = S — 2eA, and 
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where m = (mim2m3)^/^ is the geometric mean effective 
mass. 

We expand the order parameter in the harmonic oscil- 
lator states and the spherical harmonics 

oc rn—1 

A(R,k) = ^ J2 \n><n\Aim>Y{^{ek,<Pk), (4) 

n=0 m=—l 

and insert Eq (4) into the Gor'kov gap equation. The 

details are presented in the appendix. Then for B \\ 63, 
the polar and two axial PBS states are obtained from 



< n|Aio > a(f) = 0, 

(< n|Aii > ± < n|Ai _i >)a^ 
Tbn-2 <n- 2|Aii > -bn <n + 2| Ai _i >, 

where 

a(f) = [N{0)V,]-'-aiP\ 

ai") = [NiO)V2]-'-al:\ 

a(+) = [N{0)V,]- 
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6„ = 7rTV/ deUin'e d^e-^^^^'^^e-^^^^ Fn{^^2), 

(11) 

where X = p (a) corresponds to the upper (lower) row 
in Eq. (10), and 62 = eB^^vlsm^9{^), where mi2 = 
,/mim2. The bare transition temperatures are given by 
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for i = 1,2, 3, where Wd is the appropriate energy cutoff, 
and 7 1.781 is the exponential of Euler's constant. 
Defining =< n|Aii > ± < n|Ai^_i >, we obtain two 
linear equations in the two variables x^ 



2a„X„±6„_2«_2+2;„_2)+^nC^rt+2-'^n+2) = 0, (13) 

where a+ = a~ + S and S = In(^), where 
Ln{z) are the Laguerre polynomials, and Fn{z) = 



n (-2)P+V(n+2)(«+l)n! 



En 
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. With some algebraic manip- 



ulation shown in the appendix , one can show that these 
equations become a recursion relation 
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a„,n = 2a„ (1 - c„)(l - c„_2) - 5(c„ + C„_2 - 2c„C„_2), 

(17) 

and 



(18) 



It is shown in the appendix that solving this recursion 
relation equation results in a continued fraction expres- 
sion 
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By solving this equation for a given value of 6, one obtains 
hc2 for various PBS states {5 > 0, 6 < 0) the ABM state 
(5 = 0), and the CBS state (S = oo). Since the low-field 
Hc2 {T) data of Huy et al. for UCoGe suggests that it has 
uniaxial symmetry, we have restricted our study to cases 
of uniaxial anisotropy given Vi = V2 7^ V3. A full or- 
thorhombic symmetry treatment, in which Vi ^ V2 ^ V^, 
would be necessary to fit the more anomalous 5-shaped 
Hc2{T) observed by Aoki et al. at higher externally ap- 
plied fields. This will be treated in a future calculation. 



along with the inclusion of the spontaneous and field- 
induced magnetization. However, as a first approxima- 
tion, we shall neglect the effects of the magnetization, as 
in the case of URhGe, in which the magnetization affects 
Bc2(T) = ij,qHc2{T) by introducing a sharp demagneti- 
zation effect near Tc- 

In Fig. 1(a), we plotted the reduced field h^2,\\c{t) — 
2eBc2{m/mi2)vp / {2ttT^)'^ versus the reduced tempera- 
ture t = T/Tc for the polar state and for a variety of PBSs 

with —0.25 < (5 < 0, where S = ln(^^). Note that these 
PBS states all have an inherent upward curvature, but 
since > T"**, the polar state dominates near T^. How- 
ever, for —0.05 < S < there is a single kink in ft.c2,||c(^)i 
and for 6 — —0.0075, there are two kinks, due to two 
crossovers between the polar and PBS states. For 5 = 
—0.1 there is no crossover to a PBS state. In Fig. 1(b), 
we plotted hc2,±c{t) = 2eBc2{m/ y/m^^Tn^)v'j^/{2TTT^Y 
versus t = T/Tc for the CBS state and for various PBSs 
with —0.5 < (5 < 0. In this case, the CBS state domi- 
nates near to T^, but there is a crossover to the PBS state 
for —0.3 < (5 < 0, resulting in a single kink in hc2,^.c{i)■ 
Below we discuss how these kinks result in upward cur- 
vature in all three crystal axis directions, and fits to the 
UCoGe data by Huy et al. will be shown. 



FITS TO EXPERIMENTAL DATA 

Using our uniaxial anisotropic model, we fit the least 
anomalous region of the upper critical field curves, where 
the symmetry in the ab plane has not yet been bro- 
ken by the field, obtained for medium purity samples 
of UCoGe using the data of Huy et al. We have not 
included the spontaneous or field-induced magnetization 
into our calculations at this point. Future attempts to 
include the magnetization may prove to be fruitful, at 
the very least to show that it's effect on Hc2{T) is mini- 
mal, as in URhGe. In Fig. 2(a), the best fit to the data 
for the field parallel to the c-axis is shown. The graph 
shows a distinct crossover from polar (solid red curve) 
to a PBS state with & 5 = —0.0065. The experimental 
data diverges from the theoretical values for larger fields. 
The inclusion of the spontaneous magnetization might 
alter the fit, but a more plausible explanation for this 
discrepancy is a field-dependent hopping term related to 
the effective electron mass. As is discussed in Aoki et 
al., this increase in effective mass may explain the 5*- 
shaped curvature of UCoGe's upper critical field in the 
anomalous high field regime [21]. Proper consideration 
of these effects will be studied theoretically in the future. 
In Fig. 2(b), the best fit to the data for 5 = -0.185, 
which shows a distinct field-induced crossover from the 
CBS to the PBS state. 

Using extremely clean single crystal samples, Aoki et 
al. has shown that the field alignment affects the up- 
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FIG. 1: (a) Plots of K2,\\c = 2eBc2{m/mi2)vp / ver- 
sus t = T/T^ for the polar state (solid black) and for a variety 



of PBSs with 



-0.25 < ^ = in(r^"7r^'^ 

,,2 



< 



-0.001. (b) Plots 

of /ic2,xc = 2eBc2{m/^mi2m3)vl/{2-KT^f versus t = T/T^ 
for the CBS state (solid black) and for various PBSs with 
-0.5 <S< -0.1. 



ward curvature in UCoGe quite substantially, especially 
at high fields, which is likely due to the field-induced fer- 
rimagnetic state. It is therefore of great importance to 
calculate the angle dependence of the slope of the up- 
per critical field, at least close to Tc. Though the angle 
dependence of the polar state was analytically shown to 
be proportional to [cos^ 9 + sin^ 0]^-^/^ by Scharn- 
berg and Klemm [2], the angle dependence of the axial 
state has not yet been calculated analytically. This angle 
dependence may prove useful in identifying the p-wave 
superconducting state present in certain materials, and 
possibly suggest new experiments on the well known p- 
wave superconductor URhGe. 
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FIG. 2: Best fits to the data of Huy et al. for iioHc2{T) in 
medium purity UCoGe [19]. (a) i?||c. Open black diamonds: 
data. The red solid and green dashed curves are for the polar 
state and the PBS states with S = —0.065, respectively, (b) 
Data for H\\b (red crosses) and H\\a (open black circles). 
The solid black and blue dashed curves are for the CBS state 
and the PBS state with 5 = —0.185, respectively. The slopes 
at Tc were adjusted to fit the data. 



CONCLUSIONS 

We found that it is possible to fit the upward curvature 
of the Hc2 (T) data from medium-purity UCoGe using a 
crossover from the polar/CBS state to a PBS state. How- 
ever, without taking account of the spontaneous magne- 
tization, the fitting parameter 5 was different for the field 
in the c-axis direction than it was for the field in either 
the a or 6 directions. At the very least, the spontaneous 
and field-dependent magnetization should be included in 
future fits, using an anisotropic itinerant ferromagnetic 
superconductor model similar to that previously studied. 



APPENDIX 

We start with the linearized Gor'kov gap equation, Eq. 
(1), and and expand A(R, k) in the spherical harmonics 
and the harmonic oscillator basis states as in Eq. (4). 



By projecting out the Aio component one obtains 

<7i'|Aio> = 2nTN{0)V3^Y^ [ d9k'Sm9k'COS^ek' 

/"°°rf^e-2«|^„|g-ieH5^4sin'ej,,(^) 
Jo 

X < n'|e-"'"^"'^|n' >< n'|Aio >,(2G) 



where a = i£,VFy sinOk' 6^'^'°' and we have used the 
well known identity e"+^ = e'^e^e~^^°''^\ which is valid 



for 
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0. Similarly, one can find 



expressions for the other two states. 
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<n'|Aii> = 27rTA^(0)^Vi^ / dOk'sin^ek 
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It can be shown that < n'|e"" " e""|n' >= L„'(|a|^) 
and < n'|e-"'"^e"''|n' - 2 >= e-2#fe'F„,_2(|a|2) 
and < n'|e-"*"^e"''|n' -I- 2 >= 6^*'^'=' i^„' (jap) where 
Ln{x) are the Laguerre Polynomials and Fn{z) = 



5 



Ep=o p!(p+2)!(n-p)! ^0111 which we can derive the 
equations for Hc2{T). Equations (20)- (22) become 



< n'lAio > 



2TrTN{0)V3^y] f dek'smBk'Cos^ek' 
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xL„/(|a|^) <n'|Aio > 



(23) 



+ 6„-2X+_2 = a„ a;„ - 6„_2a;„_2 , (29) 
where we can let n — )• n — 2 in Eq (28) to obtain 

bn-2X^ + Q;^_2a;^_2 = bn-2X^ — Oi^_2X^_2 (30) 

a+x+ + 6„-2a;+_2 = a'x' - bn-2X~_2 (31) 

We now have two equations in the two unknowns a;+ 
and x^_2, resulting in 



< n'|Aii > 
< n'lAi _i > 



< n'|Aii > 
+ < n'|Ai _i > 



2TrTN{0)^V2Y^ J dOk'Sin^ek' 

X n'|Aii > - < n'|Ai _i 
-F„.(|a|2) <n' + 2|Ai,_i > 
+f„._2(|an < n' - 2|An >] (24) 

27rTA^(0)^Fi V / d^fc'sin^^fe. 
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ieH5^«|sin'e^,^ 
Ln'i\a\') 
X n'lAii > + < n'|Ai _i 
-F„,(|an <n' + 2|Ai,_i > 
-F„._2(|an < n' - 2|An >] (25) 



These equations may be rewritten as 



< n|Aio > ai^l2 = 0, 



(26) 



2a±a;^±6„_2(a;+_2+x„_2)+?'n(a;++2-a;„+2) = 0, (27) 

where was defined following Eq. (12) in the text. 

We note that the expression for the polar state sepa- 
rates from that for the two axial states, which comprise 
the PBS state for Vi ^¥2- In order to obtain the recur- 
sion relations for the PBS state, we begin by adding and 
subtracting both equations in Eq. (13). This leads to 
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^n-2 — ^n-2 I ^ + 



5at 



5bn-2X„ 



oit(4i-2 - bl-2 J oita!^_2 - 6,^ 



n-2 

(33) 

Letting n — >■ n + 2 in Eq. (33) and setting the two 
equations equal, we obtain 



= -2 



<5a«+2 



^ott-2 



SbnXj^_^_2 



o^n+2<^n - bl a+a+_2 - bl_2 

Sbn-2X~_2 



o^n+2<^n - bl a+a+_2 - bl_2 ' 



(34) 



where 5 = ln(Tci/Tc2), from which we obtain the recur- 
sion relation 



U — (^n,n+2X^_^2 ' ^n,n^n ~r ^n,n — 2X^-2, (35) 



where 



Sbn ,^ ^ 

an,n+2 = + — (.J- — Cn-2), 



a, 
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(36) 



(37) 



n-2 



an,n = 2a„ (1 - C„)(l - Cn-2) - S{Cn + Cn-2 - 2c„C„_2), 

(38) 

and 



(39) 



^nO.n+2 "n 



We note that there was an unfortunate typo in the ex- 
pression for Cn in Ref. (2) . 

The solution to the recursion relation can be shown to 
be a continued fraction equation involving the reduced 
field h and the reduced temperature t 
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Solutions to this equation are obtained using standard 
numerical techniques. It is important to note that the 
first order iteration differs from the exact curve by less 
than 2%. It is therefore only necessary to consider the 
continued fraction up to the 022 term, as the iteration 
thereafter converges very rapidly. 

We note that a very brief version of this work appeared 
previously [28]. 
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